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Abstract. The tetrameric Ni(II) spin cluster Ni4Moi2 has been studied by INS. The data were analyzed 
extensively in terms of a very general spin Hamiltonian, which includes antiferromagnetic Heisenberg inter- 
actions, biquadratic 2-spin and 3-spin interactions, a single-ion magnetic anisotropy, and Dzyaloshinsky- 
Moriya interactions. Some of the experimentally observed features in the INS spectra could be reproduced, 
however, one feature at 1.65 meV resisted all efforts. This supports the conclusion that the spin Hamil- 
tonian approach is not adequate to describe the magnetism in Ni4Moi2- The isotropic terms in the spin 
Hamiltonian can be obtained in a strong-coupling expansion of the Hubbard model at half-filling. Therefore 
detailed theoretical studies of the Hubbard model were undertaken, using analytical as well as numerical 
techniques. We carefully analyzed its abilities and restrictions in applications to molecular spin clusters. As 
a main result it was found that the Hubbard model is also unable to appropriately explain the magnetism 
in M4M012. Extensions of the model are also discussed. 
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1 Introduction 

In the recent years molecular nanomagnets have attracted 
huge interest because of their sometimes spectacular mag- 
netic properties. For instance, in molecules such as Mni2 
or Fes slow magnetic relaxation or even quantum tun- 
neling of the magnetization have been observed 
A general definition which includes all relevant possibili- 
ties is difficult, but in most cases molecular nanomagnets 
consist of magnetic metal ions with 3d shells and organic 
ligands. Also, in most cases the magnetism is very well 
described by assuming localized magnetic moments, such 
that metal ions with quenched orbital angular momentum 
are described simply by their total atomic spin Si (i num- 
bers the metal centers in the molecule). This yields a spin 
Hamiltonian, which typically includes a Heisenberg ex- 
change term describing the magnetic interactions between 
the metal centers in a molecule, a zero-field-splitting (ZFS) 
term describing the magnetic anisotropy due to the ligand- 
field effect, and a Zeeman term accounting for the effects 
of an applied magnetic field 0]. However, also more com- 
plicated terms may be relevant, such as Dzyaloshinsky- 
Moriya (DM) interactions, 3-spin and 4-spin exchange in- 
teractions, or higher-order ligand-ficld terms. 

Today, molecular nanomagnets comprise a large num- 
ber of magnetic molecules, and it seems fair to say that 



the larger a molecule is the more interesting it tends to 
be. Here larger means a larger number N of metal centers 
and/or larger spin lengths Si, i.e., a larger dimension of 
the Hilbert space. A larger dimension typically gives rise 
to a richer structure of the energy spectra and wave func- 
tions, and wherewith to more interesting magnetic phe- 
nomena. Also, technical aspects such as the calculation of 
magnetic observables from the spin Hamiltonian become 
quickly difficult, which establishes an interesting challenge 
by itself. 

In this context, i.e., understanding the magnetic exci- 
tations in large spin clusters, interesting molecules are for 
instance the antiferromagnetic (AFM) molecular wheels, 
in which 6 to 18 exchange-coupled metal ions form rings 

■ These wheels en- 
abled the observation of phenomena such as the rotation 
of the Neel vector [TJ]], which is related to the tower of 
states in antiferromagnetics [20,21 , quantized AFM spin 
waves 1!) 22 . or in the larger wheels quantum tunnel- 
ing of the Neel vector and the associated quantum inter- 
ference effects [23 , 24 , 25 , 26 . Numerically, for hexanuclear 
wheels the magnetism is easily calculated exactly from 
the spin Hamiltonian 27.128] . and in part also for CsFes 
with a dimension of the Hilbert space of 1.679.616 |13] , 
However, for the Feis wheel [T7], with a Hilbert space as 
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Fig. 1. (a) Ball-and-stick representation of the magnetic core 
of M4M012. Big green balls: Ni(II) ions, grey balls: Mo, small 
golden balls: O. H atoms were omitted for clarity, (b) As- 
sumed coupling paths (red/dark sticks), orientations of the 
local anisotropy axes (yellow/light arrows) and DM vectors 
(blue/dark arrows) of M4M012. 



large as ps 10 , the interpretation of the experiments re- 
quires advanced approximate techniques, which build on 
physical insight |26) . Another interesting molecule is the 
Keplerate Mo72Fe3o [25], in which 30 Fe(III) ions occupy 
the symmetric sites of an icosidodecahedron giving rise to 
strong magnetic frustration effects ,30i . Phenomena such 
as plateaus in the magnetization at 1/3 magnetization due 
to competing spin phases [3T1I32] or the presence of low- 
lying singlets [33J were observed or predicted. The Hilbert 
space is huge, of dimension rj 10 23 , hence a detailed un- 
derstanding of the magnetism in this cluster is obviously 
difficult [3J[3S]. 

In contrast to this trend to larger molecules, some- 
times even very small magnetic molecules, which at first 
sight would be discarded as trivial because of their small 
Hilbert space (which implies that "everything" can be cal- 
culated easily), may exhibit striking magnetic behavior. 
For instance, the Fe(III) dimer molecule [Fe2Fg(Et 4 N)3] 
exhibits unusual magnetization dynamics at low temper- 
atures with signatures of quantum tunneling of the mag- 
netization 36, 371135] . the unusual magnetic behavior in 
the Cu(II) tetrahedron [Cu40Cledaca4] was associated to 
phonon interactions [35] and the Ni(II) single-molecule 
magnet [Ni 4 ((hmp)(t-BuEtOH)Cl)4] allowed one to real- 
ize quantum superpositions of high-spin states [40ll411l42] . 

In this work we will study the tetrameric Ni(II) spin 
cluster [Moi 2 03o(M2-OH)ioH 2 (Ni(H 2 0)3)4], or Ni 4 Mo 12 
in short, which is another such example [43]. In M4M012 
four Ni(II) ions occupy the vertices of an almost perfect 
tetrahedron and exhibit AFM nearest-neighbor Heisen- 
berg interactions. The core of M4M012 is shown in Fig.QJa) 
The magnetism of this "simple" cluster should be unspec- 
tacular and easy to describe. Due to the AFM Heisenberg 
interactions, the ground state should belong to total spin 
S = 0, followed by a sequence of S = 1, 2, 3, 4 states with 
energies obeying the Lande rule E(S) — hAS(S + 1). In 
an applied magnetic field, this should give rise to a series 
of level crossings (LCs), where the ground state changes 



from S = to S = 1, S = 1 to S = 2, and so on. This 
sequence of LCs should in turn be detected, at low tem- 
peratures, in the magnetization curve as a sequence of 
steps at regular fields B n = nAj (g^B ) , with n = 1, 2, 3, . .. 
5J. For M4M012 magnetization steps have indeed been 
observed, but at fields of 4.5, 8.9, 20.1, and 32 T, which 
is incompatible with the Heisenberg picture |44j . This dis- 
crepancy could not be resolved by additionally introducing 
ZFS and biquadratic exchange terms in the spin Hamilto- 
nian, and a magnetic-field dependence of the exchange and 
ZFS parameters was hence proposed [33]. Subsequently, 
Kostyuchenko pointed out that 3-spin interactions should 
not be neglected [45 . Such terms can either originate 
from spin-phonon interactions as suggested in Ref. [33] 
or from electron derealization as described by a Hubbard 
model [31]. Indeed, the inclusion of such terms allowed 
Kostyuchenko to reproduce the field positions of the mag- 
netization steps. Also, starting from a Hubbard model at 
half-filling, the strengths of the Heisenberg, biquadratic, 
and 3-spin interaction terms in the spin Hamiltonian were 
obtained. According to this result, the biquadratic and 
3-spin interactions should be related by a factor of 2 (in 
our units), which was found to be in agreement with ex- 
periment as determined from a fit of the model to the 
experimental field positions. This result was interpreted 
as to indicate the superiority of the Hubbard model for 
M4M012. More recently, Klemm and Efremov presented 
an extensive analysis of the magnetization in tetrameric 
spin clusters based on a general spin Hamiltonian |46) . 
However, unfortunately, the symmetry case relevant for 
N14M012 was not considered. 

In order to better understand the unusual magnetism 
in M4M012, we undertook inelastic neutron scattering (INS) 
experiments as well as a more detailed analysis of the Hub- 
bard model, which we will present in this work. INS is 
renowned for its ability to study exchange splittings in 
magnetic clusters directly [Ill371iSJI3i[Sni[SI][i3[S3] . Our 
extensive analysis of the data in terms of a very gen- 
eral spin Hamiltonian will provide insight into the im- 
portance of the various interaction terms, but we could 
not find a parameter set which reproduces all key aspects 
of the data. This is certainly an unsatisfying outcome, 
but emphasizes the unconventionality of the magnetism 
in N14M012. It also suggests to study models going be- 
yond the spin Hamiltonian approach, such as the Hubbard 
model [54] . which allows for mobile electrons, i.e. itinerant 
magnetic moments. At half- filling the electron mobility is 
governed by the Hubbard-C/ parameter, and the model 
of localized moments is recovered in the limit of large U 
(strong-coupling limit). Our work on the Hubbard model 
is the most careful application of it to a molecular spin 
cluster to date. However, in contrast to Ref. [35], we find 
using various techniques that the Hubbard model is not 
adequate for describing the magnetism in M4M012. We 
will clarify this discrepancy, and as a byproduct resolve 
some errors in previous works. 

We like to mention that a Hubbard model descrip- 
tion of a magnetic molecule could be of interest by its 
own. For example, the low-energy spectrum might con- 
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tain additional levels not present in a pure spin model, 
which may be crucial for an appropriate theoretical de- 
scription of the magnetism in a molecule. As mentioned 
earlier, also biquadratic multiple-spin interactions such as 
2-, 3-, and 4-spin terms appear in the strong-coupling limit 
of the Hubbard model [55] . There is thus the prospect of 
identifying electron derealization as suggested by Density 
Functional Theory calculations [SS] and modeled by the 
Hubbard Hamiltonian as the physical mechanism leading 
to such interactions and giving a better understanding of 
them. 

The remainder of this article is organized as follows. 
In section [5] the experiments and a first qualitative analy- 
sis are described. Section [3] presents the phenomenological 
spin Hamiltonian that is used in section 0] to analyze the 
data with respect to various parametrizations. Section [5] 
discusses a possible description of N14M012 in terms of a 
Hubbard model. The article closes with conclusions. 



2 Experiments and Qualitative Analysis 

A fully-deuterated powder sample of M4M012 was used for 
the INS experiment. For details on synthesis and molecu- 
lar structure see [13]. The INS data were measured on 
the direct time-of-flight spectrometer IN5 at the Insti- 
tut Laue-Langevin (ILL, Grenoble, France). Spectra were 
recorded at temperatures of 2.4, 9.3, and 23 K for an in- 
cident neutron wavelength of A = 5.0 A. The energy res- 
olution at the elastic peak was 118 /ieV. The data were 
corrected for detector efficiency via a vanadium standard, 
and spectra were summed over all detector banks. 

Fig. presents the INS spectra. Four features can 
be identified on the neutron energy-loss side. At 2.4 K 
two strong peaks at about 0.4 and 0.65 meV (peaks I 
and II henceforth) are clearly visible. Their intensity de- 
creases significantly with increasing temperature. Hence, 
these peaks correspond to cold transitions and are mag- 
netic in origin. At a higher energy of about 1.7 meV (peak 
III) additional intensity is also visible at the lowest tem- 
perature. This broad feature is hence also assigned to 
a cold magnetic transition. At intermediate energies of 
about 1.1 meV (peak IV) additional intensity appears at 
higher temperatures, which we thus assign to a hot mag- 
netic peak. The sharp feature at 1.4 meV marked by an 
asterisk in Fig. [21a) is a spurion. The scattering intensity 
on the neutron energy-gain side increases with increasing 
temperature, in accord with general expectations and in 
agreement with the above assignment of the features I to 
IV, but is pretty featureless, which is consistent with the 
reduced resolution for up scattering in direct TOF instru- 
ments. We hence focus on the data on the energy-loss side 
in the following. 

In order to be able to better analyze the features, we 
subtracted a background curve from the neutron energy- 
loss data [shown as a black solid line in Fig. E^a)], which 
we generated from an educated guess. Overall the true 
contribution from the nonmagnetic scattering should be 
approximated well by our curve, but it might of course 
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Fig. 2. (a) INS spectra of N14M012 at the indicated tempera- 
tures. The black line represents the background we chose (see 
text), (b) The INS spectra with the background subtracted. 
The black lines are multi-Gauss fits with parameters given in 
Tabled] 

Table 1. Results of a multi-Gauss fit analysis of the INS spec- 
tra of Ni4Moi2 discussed in the text. The peak positions (in 
meV), line widths (FWHM, in meV), and peak intensities (in 
10~' 5 arb. units) of the fitted Gauss curves are listed. 





I 


II 


III 


IV 


energy 


0.408(3) 


0.663(4) 


1.65(1) 


1.08(5) 


FWHM 


0.15 


0.15 


0.30(2) 


0.23(2) 


intensity © 2.4 K 


2.26(4) 


1.74(4) 


1.33(5) 





intensity © 9.3 K 


0.73(2) 


0.47(2) 


0.48(3) 


0.60(3) 


intensity @ 23 K 


0.11(2) 


0.12(2) 


0.07(3) 


0.39(2) 



be incorrect regarding finer details. The background cor- 
rected spectra are plotted in Fig. E2b), showing the four 
features I to IV more clearly. Multi-Gauss fits to the peaks 
were done for each temperature, and the obtained peak 
positions, widths, and intensities are compiled in Table [1] 
The INS data indicates the energy level structure shown 
in Fig. [3] The cold peaks I, II, and III are assigned to tran- 
sitions from the ground state to three levels at 0.4, 0.65, 
and 1.7 meV. The energy of peak IV matches well the 
gap between the two lowest excited levels and the level at 
1.7 meV, which suggests to assign it accordingly. However, 
the temperature dependence of its intensity is maximal at 
around 10-15 K, such that this peak should in fact orig- 
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Fig. 3. Energy level structure as deduced from the INS data. 
The assignment of the levels to total spin quantum numbers is 
discussed in the text. 



inate from another higher lying level. The line widths of 
peak I and II are somewhat larger but close to the exper- 
imental resolution, which suggests that they are made up 
of single transitions (broadened by e.g. J strain). How- 
ever, peak III is significantly broadened, which indicates 
that a band of energy levels exists at 1.7 meV. The width 
of peak IV is also larger than the experimental resolution, 
but in view of the weakness of the signal conclusions are 
not obvious. 

A first attempt towards the interpretation of the INS 
data can be made by comparison to the magnetization 
data. In the magnetization curve at low temperatures four 
steps were observed at fields of 4.5, 8.9, 20.1, and 32 T 
[4"4"j . In a Heisenberg-exchange picture these emerge from 
the Zeeman splitting of the lowest multiplets for each to- 
tal spin S = 0, 1, 2, 3, 4, as described in the introduction. 
Since the Zeeman splitting of the involved states is known 
to be —gfisSB, the energy of the lowest spin states in zero 
field can directly be determined from the magnetization 
steps. This yields that in zero field the lowest S = 1, ...,4 
levels should be at energies of 0.6, 1.7, 4.4, and 8.5 meV 
(the ground state is S = 0). Comparing these values to 
the INS findings immediately suggests to assign the peaks 
I and II to the S = 1 level, with a ZFS of its M = and 
M — ±1 components of 0.25 meV, and peak III to the 
S = 2 level, with a weaker ZFS giving rise to an enhanced 
line width. This assignment to spin levels is also indicated 
in Fig. El 

Energy wise, this interpretation is satisfying as it is 
perfectly consistent with the observed magnetization curve. 
However, it is inconsistent with basic rules governing INS 
intensity in two instances: (1) Peak III would correspond 
to a transition from aS = 0toaS = 2 level, which is 
forbidden by the INS selection rule AS = 0, ±1. (2) The 
S = 1 levels would already be thermally populated at a 
temperature of 2.4 K (we estimate a thermal population 
of ca. 15%) and significant scattering intensity should be 
observed at this temperature, in disagreement with exper- 
iment. The above interpretation has thus to be abandoned 
as too simple. The behavior of peak III will in fact turn 
out to be a major obstacle for all the models discussed in 
this work. 



3 Phenomenological Spin Hamiltonian 

In view of the findings in previous works [44 ,45,57 and the 
above INS data, the spin Hamiltonian for M4M012 must 
include additional non-Heisenberg interactions. Here we 
consider the interactions 



T~L e x — J ^ Si ' Sjt 

i<j 
4 



4 

£ 3 = -J 3 £ (Si ■ s 3 )(s, ■ s k ), 

4 

T-Ldm = - djj ■ (Si x Sj), 
4 

Hzfs = £ Si ■ Di ■ Si, 



(la) 
(lb) 
(lc) 
(Id) 
(le) 



where H. ex describes the AFM nearest-neighbor Heisen- 
berg interactions (J < 0), ^2 and H3 are the isotropic 
fourth-order exchange terms linking 2 and 3 spins, respec- 
tively, H dm describes the antisymmetric DM interactions 
58,52|, and ~Hzfs describes the ZFS due to the local on- 
site magnetic anisotropy. To avoid confusion later on, we 
repeat that Si denotes the spin operator of the ith. Ni(II) 
ion (Si = 1 for all i). 

In principle, if going to fourth order in the isotropic ex- 
change then also a term linking 4 spins should be included. 
Its general structure is (Si ■ Sj)(Sk ■ Si) with four differ- 
ent indices i,j,k, and /. According to Ref. [45], this 4-spin 
contribution vanishes in the case of Ni4Moi2 ; we hence dis- 
regarded it in our analysis. Our discussion in section 5 will 
show that there is no reason to assume that the 4-spin in- 
teractions cancel out, but we checked that neglecting them 
does not affect our general conclusions. 

Because of the nominal high symmetry of N14M012 and 
to avoid a further increase in the number of parameters, 
we assumed that the isotropic coupling constants J, J2, 
and J3 are equal for each coupling path. The Hamilto- 
nian T~L ex + H2 + %3 can be solved analytically by using 
the vector-coupling rules of angular momenta (Kambe's 
method) [45,46^ However, the eigenvalues provided in Refs. 
and [IS] disagree. We hence performed systematic test 
calculations using exact numerical diagonalization, which 
confirmed the eigenvalues in Ref. |45j . 

According to the DM rules [S5], the DM vectors dij 
should vanish for a perfectly tetrahedral cluster. However, 
a slight distortion from this high symmetry allows for non- 
zero DM interactions, which can in fact become significant 
[591160) . Assuming a S4 symmetry, the orientations of the 
dij shown in Fig.jljb) result, with equal lengths of the DM 
vectors for all coupling paths. We therefore write dij = 
deij with appropriate unit vectors e^. 
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The on-site anisotropy terms may also be effective in a 
perfectly tetrahedral cluster, but the local anisotropy ten- 
sors are obviously severely constrained in their orientation 
and magnitude by symmetry. To avoid over-parametrization, 
we assumed an axial local anisotropy on each Ni site, each 
characterized by a tensor D = Z)diag(— |, — 5, |) in the 
local coordinate frame. The local coordinate frames are re- 
lated to the cluster coordinate frame by rotation matrices 
R{on, ipi), with the Euler angles a*, and ifi [HI]- In 
the cluster frame the local anisotropy tensors hence be- 
come Di = R(ai,$i, ipi) -D ■ R 7 \ai,$i, ipi). The local z 
axes were chosen to point radially outwards, as shown in 
Fig.[TJb). In principle, any configuration of the axes which 
is related to this one by a global rotation would also satisfy 
the tetrahedral symmetry. However, such a global rotation 
does not affect the magnetic properties of powder samples 
(we neglect the weak interplay of local anisotropy and DM 
interactions). We mention that we also considered an ad- 
ditional orthorhombic on-site anisotropy, i.e., a local ZFS 
tensor of the form D = diag(-±L> + E, -\D - E, §£>), 
but this did not lead to a significant improvement or fur- 
ther insight. We hence disregarded this term. 

In an applied magnetic field, the Zeeman term 

4 

U B = VbqB ■ Si (2) 

i=l 

is additionally present, where g is typically on the order 
of 2.3 for Ni(II) ions. We have neglected an anisotropy of 
the g tensor, as it is expected to be small and irrelevant. 

For a given spin Hamiltonian, the magnetization curve 
as well as the INS spectra were calculated numerically 
from the eigenpairs obtained from a full exact numerical 
diagonalization. For the magnetization curves, the powder 
average was done by numerically averaging over a grid of 
magnetic field orientations. The powder INS spectra were 
calculated using the formulae given in Refs. [62ll63j . 




energy (meV) 

Fig. 4. Experimental (dots) and simulated (lines) INS spectra 
for H = U ex + Uzfs with J = -6.5 K and D = 3.25 K. 

Table 2. Parameters of the three models discussed in Ref. |44] , 



model 


J a 


J b 


JS 


4 


D 


a 


-6.4 K 


-6.4 K 


3.2 K 


K 


-1.0 K 


b 


-6.4 K 


-6.2 K 


3.0 K 


K 


-3.2 K 


c 


-8.4 K 


-6.2 K 


0.32 K 


0.78 K 


-8.9 K 



terms in Ni4Moi2- In order to be reasonably consistent 
with the experimental data (magnetization and INS), we 
observed that D cannot be very large, such that the mix- 
ing effect is not very strong. Accordingly, the simulated 
intensity for peak III is way too weak. The best simula- 
tion of the INS spectra was obtained for J = —6.5 K and 
D = 3.25 K, see Fig. |U The peak positions are well re- 
produced, but the discrepancies as regards the scattering 
intensity are obvious. Most noteworthy, the intensity of 
peak III is strongly underestimated and exhibits a totally 
wrong temperature dependence. We couldn't find any pa- 
rameter set for which these discrepancies did not arise. 



4 Analysis 

4.1 Heisenberg exchange plus on-site magnetic 
anisotropy 

The qualitative interpretation of the INS spectra in terms 
of a dominant Heisenberg interaction has shown that peak 
III would violate the INS selection rule AS = 0, ±1. How- 
ever, for Ni(II) ions it is well known that they may exhibit 
on-site anisotropies as large as several ten K [64]. Hence, 
the possibility arises that in Ni4Moi2 the anisotropy con- 
stant D is on the order of or even larger than the exchange 
coupling J, which would give rise to strong mixing of spin 
levels, such that S would cease to be a good quantum num- 
ber [63,65,66 . Then also AS = 0,±1 would cease to be 
a good selection rule. In a hope that this mechanism may 
explain the significant INS intensity of peak III, we did ex- 
tensive simulations for the Hamiltonian H — Wex+'HzFS, 
scanning the whole parameter regime for J and D. This 
model accounts for the two most significant interaction 



4.2 The models of Schnack et al. 



Based on the magnetic susceptibility and magnetiza- 
tion curves, three models were suggested for M4M012 in 
Ref. 44 . In these models additionally a structural dis- 
tortion was accounted for via different coupling constants 
for the isotropic exchange, J12 = J2.3 = ^13 = J a and 
J14 = J24 = J34 = J b ' 1 and similarly for the 2-spin inter- 
actions, introducing parameters and J\. The parame- 
ters of the three models are reproduced in Table [5] (in our 
units), where the suggested magnetic field dependence of 
some of the parameters is irrelevant as our INS experi- 
ment was done in zero field. The simulated INS spectra 
are compared to the experimental data in Fig. [5] Unfor- 
tunately, none of these models can reproduce the experi- 
mental INS spectra. The effects of DM interactions and/or 
tilted anisotropy tensors should be even weaker (the effect 
of DM interactions will be discussed in more detail in sec- 
tion 4.4). Hence, a deviation from the assumed tetrahedral 
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Fig. 5. Experimental (dots) and simulated (lines) INS spectra 
for the three models suggested in Ref. [44], see also text and 
Tabled 
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Fig. 6. Experimental (dots) and simulated (lines) INS spectra 
for the model suggested in Ref. [IS], TL = H ex + H2 + H3 with 
J = -8.82 K, J 2 = -1.08 K, and J 3 = -0.53 K. 



obtained relations are erroneous. Nonetheless, used as a 
phenomenological model, H = H ex + H2 + H3 + Hb very 
successfully describes the magnetization data for the pa- 
rameters J = -8.82 K, J 2 = -1.08 K, and J 3 = -0.53 K 
which were obtained from a fit to the experimentally de- 
termined crossing fields [45) . However, this model does not 
reproduce the experimental INS spectra well, see Fig. [B] 
The peak positions may be considered to be acceptable 
in view of the simplicity of the model, but we again ob- 
serve the same discrepancy as before, i.e., the intensity of 
peak III is strongly underestimated and exhibits a wrong 
temperature dependence. 

Because ZFS can be significant for Ni(II) ions, we ex- 
tended the model to % = ii ex + H 2 + ^3 + ~Hzf s +Hb 
and ran a fit to the data, in which we used the magnetiza- 
tion and INS data simultaneously. The resulting best-fit 
curves are shown in Fig. [7] The magnetization is well re- 
produced, but the fit to the INS data shows again the 
"peak III discrepancy". Moreover, now also peak IV is 
incorrectly reproduced both as regards its intensity and 
temperature dependence. Apparently, also this extended 
model cannot satisfactorily account for peak III. 



4.4 Systematic scan of the parameter regime 



symmetry is apparently not at the heart of the unusual 
magnetism in Ni4Moi2- 

4.3 The model of Kostyuchenko 

Kostyuchenko showed recently that the positions of the 
steps in the magnetization curve can be reproduced very 

well by the model % = Hex + % + H3 +~Hb, which in 
particular includes the 3-spin interactions [15] . He also de- 
rived this spin model from the strong-coupling limit of a 
Hubbard model at half-filling, and found that the cou- 
pling strengths J, J2, and J3 are not independent. Un- 
fortunately, as our analysis in section 5 will reveal, the 



In a last effort to explain the magnetism in M4M012 in 
terms of a spin Hamiltonian we systematically considered 
the Hamiltonian H = Hex + H2 + H3 + Hzfs + Hdm, 
which includes all terms considered so far [see Eqs. (JTJl]. 
We started from simulations of the INS spectra for the 
Heisenberg Hamiltonian, Eq. (fTa)) , and then extended the 
model by including the other terms step by step. In every 
step, the influence of the newly added term was analyzed. 
We will not discuss all details of our findings here, but 
just mention some key observations. 

With the Heisenberg term, Eq. (fTa|) . alone, no agree- 
ment between data and simulation could be obtained. We 
hence chose J = —8.82 K. The effects of the further terms 
on the INS spectra are presented in Fig. [5] Each panel 
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J 2 (K) J 3 (K) D/J d/J 

Fig. 8. Simulated INS intensity as function of energy and a selected magnetic parameter. The strength of the scattering 
intensity is represented by the color (red = strong, blue = weak), (a) Intensity vs. energy and J2 at 2.4 K. (b) Intensity vs. 
energy and J3 at 2.4 K (J2 = — 1.08K). (c) Intensity vs. energy and D at 2.4 K (J3 = — 0.53K). (d) Intensity vs. energy and d 
at 9.3 K (D — 3.5 K). As a guide to the eyes, the black lines indicate the energies of the measured peaks I, II, III, and IV. 




energy (meV) 

Fig. 7. (a) Experimental (black dots) and simulated (green 
line) magnetization curve and (b) experimental (dots) and sim- 
ulated (lines) INS spectra for "H = Hcx+H2 + 'H3 + 'Hzfs + 'Hb 
with J = -9.45 K, J 2 = -1.41 K, J 3 = -0.66 K, D = 6.26 K, 
and g = 2.35. 



presents the INS intensity at one temperature, plotted 
as function of energy and the magnetic parameter under 
consideration, with the scattering intensity represented by 
color. The dependence of the INS intensity at 2.4 K on the 
strength of the 2-spin interaction J2 is shown in Fig.[5Ja). 



Two pronounced peaks at ca. 0.5 and 0.8 meV and a weak 
peak at ca. 1.8 meV are visible. They are close to the ex- 
perimental energies of features I, II, and III, but exhibit a 
rather weak dependence on J2 . Hence, the 2-spin interac- 
tion is important, but its exact strength is not constrained 
much. We chose J2 = —1.08 K. The simulations with ad- 
ditional 3-spin interactions are shown in Fig. |8]Jb). The 
3-spin interactions have a pretty strong effect on the INS 
spectrum. The two low-energy features, which were visible 
in Fig.[8ja), are now highly entangled, and the high-energy 
feature becomes very weak. The 3-spin interactions were 
crucial for explaining the magnetization data, but as re- 
gards INS they actually have a counter productive effect. 
If it were not for the magnetization, one would rather 
abandon them in an interpretation of the INS data. We 
chose J3 = —0.53 K. Next the ZFS term was included; the 
simulations for varying D values are shown in Fig. ISJc). 
As expected, the ZFS term has a significant effect, and 
for D values in between — 0.5J and — 0.1J both peaks I 
and II are well reproduced. However, as discussed before, 
even for large D values the mixing of spin levels is not 
strong enough to produce a significant INS intensity at the 
higher energies, i.e., peak III is not reproduced. Further- 
more, at 2.4 K there is significant intensity in the energy 
region 0.8 to 1.2 meV, in contrast to experiment. We chose 
D = -0.4 J = 3.5 K. Finally, the effect of the DM interac- 
tions on the INS spectrum at 9.3 K is shown in Fig. [Sfd). 
The DM interactions could in principle also be responsi- 
ble for a mixing of spin levels and hence a violation of the 
INS selection rule, but apparently they have no essential 
effect. We therefore conclude that DM interactions are not 
relevant for explaining the key unexplained aspects of the 
magnetism in N14M012. 

The best simulation obtained in this approach is com- 
pared to the experiment in Fig. [HI Peaks I and II are repro- 
duced reasonably well. Also the temperature dependence 
of peak IV is by and large correctly obtained (increasing 
intensity with increasing temperature). However, as for all 
models discussed in this work, the intensity and tempera- 
ture dependence of peak III is not reproduced. 
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energy (meV) 

Fig. 9. Experimental (dots) and simulated (lines) INS spectra 
for our best-fit model, H — H ex + H2 + Ha + Hzfs with 
J = -8.82 K, J 2 = -1.08 K, J 3 = -0.56 K, and D = 3.5 K. 

To conclude this section, all our efforts to reproduce 
the experimental INS data with a spin Hamiltonian were 
unsuccessful in the sense that key aspects of the data, 
such as the magnitude of the intensity of peak III and its 
temperature dependence, could not be satisfactorily re- 
produced. The analysis, however, provided some insight 
into which terms should be most important. For example, 
deviations from the tetrahedral symmetry and DM inter- 
actions appear to be irrelevant. Considering the extensive 
yet unsuccessful spin-Hamiltonian-bascd efforts presented 
here as well as in Ref. [H], one may wonder about the 
apparent inability of the spin Hamiltonian approach to 
account for the magnetism in N14M012. 



5 Hubbard Model Description 

The difficulties in interpreting the experimental data by a 
spin Hamiltonian suggest to reanalyze its basis. A possible 
approach is to start with the more fundamental Hubbard 
model at half-filling, from which spin exchange interac- 
tions are obtained with a standard argument in a large 
U expansion. In fact, a Hubbard model was recently pro- 
posed for M4M012 [45]. However, as we will show in this 
section, this model is not adequate for describing the low- 
tempcrature thermodynamic and spectroscopic properties 
of the molecule either. 

As presented in Ref. [IS], the motivation for applying 
a Hubbard model in the first place comes from the octahe- 
dral surroundings of the Ni ions which cause a splitting of 
their 3d levels into three lower-lying t-i g and two higher- 
lying e g orbitals. According to Hund's rules, the t 2g or- 
bitals are fully occupied, whereas the e g orbitals are only 
singly occupied. The proposed Hubbard model is then for- 
mulated using a standard (one-band) Hamiltonian |54) . 

L L 

Tin = tgpcl^cp^ + U}, n at fh a> i + gpL B S ■ B , (3) 

afS.a a=l 




Fig. 10. Schematic representation of the lattice and the hop- 
ping paths for the Hubbard model Hh (cf. Fig. 2 of Ref. [45] 1. 
The small, numbered spheres correspond to Hubbard lattice 
sites, which are assigned to the four Ni ions as indicated by 
the large, green circles. Possible jumps of electrons are visual- 
ized by lines: inter-ion hoppings (dashed/red) occur with equal 
strength t and intra- ion hoppings (dotted/blue) with strength 
t a , but the latter are excluded and set to zero. 

with parameters N e = L = 8 (half-filling), where L is the 
number of lattice sites and N e the number of electrons. 
The c aa (cu^cr) are the fermionic creation (destruction) 
operators for an electron at lattice site a with spin projec- 
tion a =t>4o an d n a , a = c a „c a ^. The local spin density 
at site a can be written as s a = '}2, IJT c a a a aT c a ^ T , with 
the vector of Pauli matrices a. The total spin operator 
then reads S — J2 a The possible jumps of electrons 
between lattice sites are characterized by the hopping pa- 
rameters t a p, as illustrated in Fig. [TU] In the considered 
model all hopping matrix elements corresponding to in- 
ter-ion jumps are equal to t, and those corresponding to 
mfra-ion jumps are equal to t a . The parameter t a is set 
to zero, as in Ref. [45] . 

In this section we will focus on the steps in the mag- 
netization curve as predicted by the Hamiltonian Hh- 
At this point we should stress that in Ref. [45] not the 
Hubbard model itself was analyzed, but the spin model 
Hex +U2 + H 3 + Ub [Eqs. (TliJ-([lc]), and ©], which is 
supposed to approximate the full Hamiltonian Hh - In con- 
trast, here we directly study the Hubbard model by using 
numerical exact diagonalization to calculate the eigenval- 
ues of Hh- The parameters t and U were then fitted to 
the experimentally determined crossing fields by employ- 
ing a numerical minimization routine without imposing 
any constraints on the parameter values. Temperature was 
set to T — 0, and the g factor to g = 2.25, as in Ref. [45 . 
Despite testing a lot of initial values and initial search 
directions for t and U, we only found a single, unsatisfac- 
tory fit which obeys U > 0, namely t — 3.91 meV and 
U = 64.6 meV. The predicted low-temperature magne- 
tization curve is shown in Fig. [TT] The x 2 value of the 
fit to the crossing fields is about 6.25 T 2 , which is ap- 
proximately two orders of magnitude worse than the re- 
sult obtained with the spin model Ti ex + y.2 + H3 + T~Lb 
[45] . Considering the stability of our fit result regardless 
of the starting conditions, and taking into account that 
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Fig. 11. Magnetization of the Hubbard model Hh as a func- 
tion of the applied magnetic field B at T = 0.01K for the 
parameters t — 3.91 meV, U = 64.6 meV, and g = 2.25 
(red/dark). The steps occur at magnetic fields of 5.0, 10.7, 
18.5, and 32.3 T. For comparison, also the fit result of Ref. [35] 
is plotted (green/light), which was obtained with the Hamilto- 
nian H ex + Hi + Hz + Hb (see section |4~3)) . The positions of 
the experimentally determined crossing fields are indicated by 
vertical lines. 



the parameter space of the Hubbard model Hh is only 
two-dimensional, we have to conclude that Hh cannot 
explain the low-temperature magnetization in N14M012. 
Accordingly, it does not represent an adequate model for 
M4M012, and thus does not capture the physical origin 
of the non-Heisenberg interactions in the molecule. This 
clearly contradicts one of the main findings of Ref. [35]. 
In the following we will hence substantiate our conclusion 
and clarify the situation. As a byproduct, we will iden- 
tify an error in the analytical derivation of the spin model 
Hex + H2+H3+H.B in Ref. [45]. 

Before doing so, we should first comment on the pro- 
posed Hubbard model Hh and point out some of its short- 
comings and possible modifications. As far as we under- 
stand it, the intra- ion hopping terms were neglected in 
Ref. |45) on the grounds that the Pauli principle prohibits 
a hopping between lattice sites which are occupied by elec- 
trons with equal spin projection. However, this argument 
is valid only in a classical picture of a spin-1 state. Quan- 
tum mechanically there is no a priori reason for neglecting 
the intra-ion hopping, because a level with s = 1 can also 
have a magnetic quantum number of m = 0. Neglecting 
the intra-ion hopping is hence an inadequate method to 
incorporate the known ground-state spin configuration of 
the Ni ions into the Hubbard model (see below). We have 
checked for a number of parameter sets that a non-zero 
value of t a can influence the low-temperature magneti- 
zation curve by shifting the position of the ground state 
energies in subspaces with fixed magnetic quantum num- 
ber. In our fits to the magnetization data we therefore also 
added t a as a third independent fit parameter. However, 
this did not lead to any improvement. Furthermore, ac- 
cording to the motivation for using a Hubbard model as 
presented before, one should actually use a two-band Hub- 
bard model because every Ni ion provides two magnetic 



orbitals. In the usual derivation of the Hubbard model (see 
e.g. Ref. [67]) this leads to a variety of additional inter- 
action terms (intra- ion Coulomb repulsion, mira-ion ex- 
change interaction, pair jumps, and correlated hopping). 
Of these we considered only the terms which are usually 
studied in the literature, i.e., longer-range Coulomb repul- 
sion and Heisenberg exchange (extended Hubbard model). 
The Heisenberg exchange term provides a proper means 
to handle the spin state of the Ni ions, through the use of 
a ferromagnetic intra- ion coupling (Hund's rule coupling). 
Perhaps not surprisingly, the inclusion of further interac- 
tion terms allowed for a better fit of the magnetization 
data and improved the x 2 value by about one order of 
magnitude. However, this is still one order of magnitude 
worse than the fit result obtained in Ref. [IS]. We thus 
do not believe that such generalizations of the Hubbard 
model Hh lead in the right direction. 

In order to resolve the apparent contradiction between 
our conclusion and that in Ref. |45] . we reanalyzed the 
strong-coupling limit U 3> \t\ of the Hamiltonian Hh- 
In this limit the hopping term can be treated as a per- 
turbation, which at half-filling leads to an effective spin 
model H s whose energy eigenvalues are supposed to ap- 
proximate the low-energy spectrum of the full Hubbard 
model [6? . It should be noted that the effective Hamil- 
tonian H s is a spin-1/2 model, consisting of 8 spins in 
the case of N14M012. We calculated Hs analytically up to 
order 0(t 4 U~ 3 ) using two approaches. Starting from the 
results of the so called Canonical Transformation [55,68, 
Hubbard-X operators [69 were employed to rewrite the ef- 
fective Hamiltonian in the form of a pure spin model, as in 
Ref. [70J . Alternatively, one can use the general results of 
Takahashi [7T] and MacDonald et al. [75] to directly find 
the effective spin Hamiltonian, which considerably speeds 
up the whole calculation. Both approaches lead to equiv- 
alent results. Denoting the spin Hamiltonian up to order 

Oit 2 !/- 1 ) as Hi 2) , and the one up to order 0{t 4 lJ- 3 ) as 
H s 4 \ we obtained the following expressions: 



H^ 



«i 4 ' = 



-j (2) S « ' S - — + 9^bS ■ B , (4a) 

a</3 

8 4 

-j (4) ^2s a - Sfs - j ion ^ ««, • S /3i ( 4b ) 
a</3 i=l 

-fa [^ a ' ^ ' Sff ) 

(a/97*) 

+ (s a ■ Sg) (sp ■ S 7 ) - (s a ■ S 7 ) (sp ■ Ss) 

(2At 2 78t 4s 



\ U U 3 



9HbS ■ B . 



The terms 



-fa 



Ea</3 ■ h P and "J (4) ELfl 5 



de- 



scribe the inter-ion interactions, where the sites a and /3 
belong to different ions; the term —ji on Yl%=i ' de- 
scribes the intra-ion interactions, where the sites ct^ and 
Pi belong to the same Ni ion i [(c^,/^) = (1,2), (3,4), 
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(5,6), (7,8) for i — 1,2,3,4]. (a^jS) denotes the tetragon 
with vertices a, /?, 7, and 5, which have to be connected 
by hopping paths such that a is connected to /3, ft to 7, 
7 to 8, and S to a (see also Fig. ITU1) . The Zeeman term 
gfigS ■ B is unaffected in the perturbation theory as the 
total spin S commutes with the hopping term. The cou- 
pling parameters were obtained as 



,■(2) _ 



(4) 



4*2 
U ' 

At 2 92t 4 

W + W' 



34 



"773" 



36t 4 
773" 



(5a) 
(5b) 



Corrections to the spin Hamiltonian T-L^ and the given 
parameters are of order 0(t 6 U~ 5 ) [71]. Even though intra- 
ion hoppings are not included in the model, a non-zero 
and ferromagnetic coupling between the spins belonging 
to the same Ni ion is obtained, ji on > 0. However, it is 
of order 0(t 4 U~ 3 ) and thus usually weak. We mention 
that the energy of the ferromagnetic ground state with 
\M\ = \M max \ is zero for both Hh and its effective spin 
Hamiltonians [the corresponding constants were explicitly 
retained in Eqs. {3}]. 

In order to evaluate the accuracy of our analytical re- 
sults for Hs and Hi, we calculated the energy eigen- 
values of the different models for a number of parameter 
sets t and U via numerical exact diagonalization and com- 
pared the obtained energy spectra. For the ratio U/t = 20 
the energy spectra are presented in Fig. [T2] (t = 1). Ap- 

parently, the spectrum of H s clearly deviates from the 
exact spectrum, in particular at low energies, whereas the 
fourth-order Hamiltonian H^ gives a much better ap- 
proximation. We found that for U/t = 30 the low-energy 
spectrum of the Hubbard model Hh is already very well 
approximated by the eigenvalues of H s ■ However, we also 
observed that for a ratio U/t = 10, which was inferred 
and identified with the strong-coupling limit U 3> |i| in 
Ref. [3S], the exact and fourth-order spectrum have little 
in common. A link between exact and approximative lev- 
els may be established starting with U/t « 15, but at this 
ratio clear discrepancies are still present, especially with 
respect to the ground state energy. In the following we will 
only consider H^ . 

Having established an effective spin Hamiltonian for 
the electron spins s Q , we now discuss how to derive an 
effective spin model in the space of the Ni spins Si. The 
total spin of the Ni ion i is Si = s ai +8p t . Since the oper- 
ators Si commute with Hs , the quantum numbers Si are 
good quantum numbers. The part of the energy spectrum 
with Si = 1 for all i can thus be extracted by adding a 
ferromagnetic coupling between the electron spins s ai and 

s / 3 i residing on the same Ni ion to Hs (Hund's rule cou- 
pling). Technically, this can be done by setting j ion — > 00 
(energy levels with Si = for at least one i are lifted 
to arbitrarily high energies). However, it is important to 
note that this approach is only possible within the effec- 
tive spin model. In the case of the Hubbard model Hh 
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Fig. 12. Numerical comparison of the low-energy spectrum 
of the Hubbard model Hh (red/dark, right) and the complete 
spectrum of the spin model H^ (blue/dark, left) and Hi 4 - 1 
(green/light, middle) for the parameters t = 1 and U = 20. 
The energy eigenvalues are plotted vs. the absolute value of 
the 'total magnetic quantum number \M\. The product state 
with \M\ — 4 is an eigenstate of all three models with energy 
E = 0. 



the Si are not good quantum numbers since any opera- 
tor which is specific to particular Hubbard lattice sites, 
like Si = s ai + Sfc, cannot be a conserved quantity in an 
itinerant model. Adding a Hund's rule coupling directly 
to T-Lh would thus make the configurations with Si = 1 
energetically favorable, but it would also change the struc- 
ture of the whole energy spectrum. This means that it is 
incorrect to state that the spin model which is obtained 

by adding such a term to Tii^ follows from the original 
Hubbard model. 

To carry out the limit ji on — > 00, the operators Si can 
be interpreted as pure spin-1 operators in a low-energy 
theory. The resulting spin model then becomes (using no- 
tation introduced in section 3) 



l~i — l~ic 



H2+H3 + H! 



(6a) 



— J4 (Si ■ Sj)(Sk ■ Si) — 

(ijki)er 

with the coupling parameters 



/24i 2 



V u 



288t 4 



u 



\mt 4 



J 3 = ^ 



40t 4 
IP 



(6b) 



and r = {(1,2, 3, 4), (1,3, 2, 4), (1,4, 2, 3)}. We men- 
tion that introducing an mira-ion hopping of strength t a 
= t (instead of t a — 0) in "Hh does not change the Hamilto- 
nian Eq. ([5]) (but certainly changes and H^), which 
is due to the high symmetry of the present case. These 
results finally clarify our different findings for the magne- 
tization steps when calculated with the Hubbard model 
%h and the spin model % ex + H2 + H3 +Hb, and should 
be compared to those in Ref. [45]. One notable difference 
is that the Hamiltonian Eq. © comprises 4-spin inter- 
actions. Starting from Eq. (|4bl) . it can easily be verified 
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that such terms have to appear in the projection of Hs 
onto the subspace with Si = 1 for all i. Furthermore, we 
find that the parameters Ji and J3 are of equal strength, 
whereas in Ref. [45] they differ by a factor of 2. Because the 
analysis of the experimental magnetization data with the 
spin model Hex + H2 + Hz + Hb yielded such a factor of 2 
and the inferred ratio U/t — 10, it was concluded that the 
Hubbard model provides a valid description of the mag- 
netism in N14M012 [73]. However, our results demonstrate 
that such a conclusion is not supported. 



6 Conclusion 

In conclusion, we presented inelastic neutron scattering 
data for the M4M012 molecule. An extensive analysis of 
these data in terms of a phenomenological spin Hamilto- 
nian did not lead to a satisfactory description of its mag- 
netism. A similar observation was made in previous works, 
mainly based on magnetization data |44] . Confirming this 
observation by complementary spectroscopic data, in our 
opinion, significantly furthers the idea that for N14M012 
the spin Hamiltonian approach is indeed inadequate. As 
an alternative and more fundamental model the Hubbard 
model comes to mind, which we hence analyzed in great 
detail. We have shown that the first-guess Hubbard model 
Hh cannot explain the magnetization data of M4M012 as 
it predicts incorrect crossing fields. Obvious extensions of 
the model did not resolve the issue satisfactorily either. 
We then studied the strong-coupling limit of Hh in or- 
der to better understand these results. According to our 
calculations, the observation in Ref. [?S] that the mag- 
netization curve can be nicely fitted with the spin model 
H ex + Hi + H3 + Hb does not imply that the Hubbard 
model is a suitable microscopic model for N14M012. In this 
work we have in fact shown that the opposite is true. 
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